function f = z_high2(zh,beta,k,R,q_bar)
% Upper bound of branches leading to bunching: proportional
f = exp(zh + beta*R - 1)/beta*(1 - k) -(R - (q_bar - zh)/beta)*exp(q_bar); 
end